Stick-Slip Sliding of Water Drops on Chemically Heterogeneous Surfaces 
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We present a comprehensive study of water drops sliding down chemically heterogeneous surfaces 
formed by a periodic pattern of alternating hydrophobic and hydrophilic stripes. Drops are found 
to undergo a stick-slip motion whose average speed is an order of magnitude smaller than that 
measured on a homogeneous surface having the same static contact angle. This motion is the result 
of the periodic deformations of the drop interface when crossing the stripes. Numerical simulations 
confirm this view and are used to elucidate the principles underlying the experimental observations. 
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Controlling and predicting the equilibrium and dynam- 
ical properties of drops on chemically heterogeneous sur- 
faces is a major scientific challenge, relevant not only 
for fundamental research [ll-Q but also for an ample va- 
riety of applications, particularly in droplet based mi- 
crofluidics, a rapidly growing interdisciplinary field of re- 
search combining soft matter physics, biochemistry and 
microsystems engineering [5]. Understanding the liquid 
dynamics is essential to design new smart coatings which, 
for example .guide the wetting drops along certain di- 
rections p, |7j. When drops are deposited on a sur- 
face functionalized with stripes of alternating wettability, 
they may assume elongated shapes, which are character- 
ized by different contact angles in the directions perpen- 
dicular and parallel to the stripes. This morphological 
anisotropy has been the object of intense scrutiny in a 
variety of situations (il-[l2|. Studies about the sliding of 
drops on striped surfaces also report an anisotropic dy- 
namics: drops move faster along the alternating stripes 
than across them 0, E3 and velocity fluctuations take 
place a long th is latter direction [T^j. Numerical simu- 



lations 14 - 16| and theoretical studies 17-19| also point 
to the occurrence of a stick-slip motion when the drop 
crosses stripes of different wettability. The observed 
anisotropics are usually explained in terms of the energy 
barrier the drop must overcome in the direction perpen- 
dicular to the stripes, while it is ideally free to move in 
the direction parallel to them. 

Hereafter, we present the first experimental evidence of 
a stick-slip motion of water drops sliding down heteroge- 
neous hydrophilic/hydrophobic surfaces having parallel 
stripes with a large wettability contrast of about 70°. We 
find that the mean speed of the sliding drop is strongly 
affected by the patterning details, with a slowing down 
that can reach up to a factor 10 with respect to the corre- 
sponding homogeneous coating with the same static con- 
tact angle. We also show that such a big change in the 
dynamical evolution is due to a different balancing be- 
tween capillary, viscous and body forces in the presence 
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FIG. 1: (Color online) a) Condensation figure on a hetero- 
geneous surface: big (small) water drops form on hydrophilic 
glass (hydrophobic OTS). Parallel stripes with periodicity W 
can be easily distinguished according to the size and distri- 
bution of the drops; b) Static profile of a water drop on the 
heterogeneous surface; c) Lateral view (top right) and con- 
tact line view (bottom left) of the sliding drops on inclined 
substrates; d) Typical Stick-Slip time law for the front and 
rear contact points of the sliding drop. 



of stick-slip. 

Chemically heterogeneous surfaces are produced 
through micro contact printing starting from the realiza- 
tion of masters with rectangular grooves by photolithog- 
raphy. Such masters are replicated in polydimethilsilox- 
ane (PDMS) to obtain the stamp [20| for the printing 
of a solution of octadecyltrichlorosilane (OTS) in toluene 
on a glass substrate. The result is a surface presenting 
alternating hydrophobic (OTS regions) and hydrophilic 
stripes (uncoated glass regions). Sample characteriza- 
tion is performed through atomic force microscopy and 
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the condensation figure method as shown in Fig. [IJl. Par- 
allel stripes of equal width and different wettability can 
be clearly evinced having a surface roughness of less than 
10 nm measured over an area of 10 /am x 10 /am. The 
resulting pattern has a periodicity W ~ 200 /im. The 
printed pattern is also analyzed in terms of contact angle 
measurements through the Cassie-Baxter equation [1[. 
Figure [TJd shows the profile of a drop deposited on the 
heterogeneous surface. Due to the size of the stripes and 
the drop volume, the effects of morphological anisotropics 
1(1 12| are not very pronounced, yielding an equilibrium 
contact angle parallel to the stripes #het,|| = 72° ± 3°, 
which fits well with the Cassie prediction #h e t = 75° ± 2° 
calculated from the pattern details. 

Water drops of the desired volume V (~ 30 fA) are 
gently deposited on the sample surface through a ver- 
tical syringe pump. The sample is placed on a tilt able 
support, whose inclination can be set with a 0.1° accu- 
racy. A mirror mounted under the sample stage at 45° 
with respect to the surface allows to simultaneously view 
the contact line and the lateral side of the drop [2l|, as 
shown in Fig.[Tfc. The drop is lightened by two LED back 
lights: one of them is placed in front of the camera and 
the other is fixed above the sample stage, to illuminate 
the drop from the top. Experimentally, we first fix the 
inclination of the plane and then deposit the drop. Ac- 
quired images, where drops appear dark on a light back- 
ground, are analyzed through a custom made program 

Fig. [T]i shows the positions of the front and rear contact 
points of a drop sliding down the heterogeneous surface 
perpendicularly to the stripes. This is more clearly seen 
in the top panel of Fig. [2j which shows the drop dynamics 
in a time window expanded by a factor ^10 with respect 
to the time scale of Fig. [T]i. The drop clearly advances 
with a stick-slip behavior, with jumps of the order of 
the pattern periodicity W. At the beginning of the time 
frame considered, the rear contact line is pinned while 
the front line slowly advances until around t ~ 0.25 s 
when it suddenly jumps forward. In correspondence, a 
drop in 9 a is observed (see the middle panel of Fig. [2]). 
Subsequently, the front contact line hardly moves, while 
the rear line at t w 0.4 s jumps by a distance equal to W. 
The process then repeats itself. This peculiar dynam- 
ics has never been observed experimentally. In previous 
studies this stick-slip motion has not been seen proba- 
bly because the sliding occurred across stripes having a 
small wettability contrast [13| . A stick-slip has also been 
reported in the spreading of a drop of increasing volume 
along the direction orthogonal to the stripes while in the 
other direction the contact line moves at a constant speed 

To better understand this motion, we have repeated 
the measurement with a glass surface coated with a ho- 
mogeneous OTS layer deposited from the vapor phase 
which presents a static contact angle # eq ~ 71° ±2°, very 
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FIG. 2: (Color online). Top: time evolution of the position of 
the rear and the front contact points of a drop sliding down 
the heterogeneous surface tilted by a=56°. Snapshots show 
the different contact angles before and after the jumps. Mid- 
dle: variation of the advancing contact angle 6a of the same 
drop and of a different one sliding down an homogeneous sur- 
face with a #eq=71° inclined by the same a. Bottom: corre- 
sponding evolution of the frontal and rear contact points on 
heterogeneous and homogeneous surfaces. 



close to the Cassie angle of the heterogeneous surface. 
Furthermore, the inclination angle of the OTS surface 
and the water drop volume are the same as in the case of 
the heterogeneous surface. Statically, then, the two sys- 
tems look pretty much the same, but dynamically their 
behavior is quite different. On the homogeneous sub- 
strate, as expected, the motion is uniform as shown in 
the bottom panel of Fig. [2] and 9 a practically does not 
change (see middle panel). Moreover, the drop veloc- 
ity on the homogeneous surface is about 10 times larger 
than the average velocity on the heterogeneous surface 
(see bottom panel). 

This physical picture is corroborated with the use of 
numerical simulations based on the lattice Boltzmann 
models [UHsl, which allow to simulate the continuity 
and the Navier-Stokes equations with capillary forces at 
the interfaces. A detailed description of the model can be 

We shall investigate 



found in other publications [27H2S 
a two-dimensional situation: this is the optimal choice 
to retain the essential features of the stick-slip dynamics 
Hi 15, ll| and, at the same time, achieve a reasonably 



large resolution to resolve the hydro dynamic equations 
inside the drop. At variance with respect to other nu- 
merical simulations done in the quasistatic limit 14, li| 
or referring to single contact line motion (l6j . we will 
show explicitly how the different terms in the equations 
of motion balance in the stick-slip dynamics. The output 
of the numerical simulations (i.e. density, velocity, forces) 
is integrated over the drop volume and made dimension- 
less with respect to the time-length scales of the problem. 
We reproduce the following force-balance equation 



Ma(t) = D(t)+F g + F ca , I >(t) 



(1) 



where a(t) is the acceleration of the drop of mass M and 
F g is the gravitational force. The term F cap accounts for 
the internal pressure forces due to the curvature distor- 
tion of the interface (F cap = for a spherical interface) 
while the function D(t) quantifies the dissipation inside 
the drop. The simulation parameters are chosen in such 
a way to reproduce the experimental conditions investi- 
gated in Fig. El The observed stick-slip motion is cap- 
tured in the numerical simulations as shown in the top 
graph of Fig. El which displays the time evolution of the 
positions of the front and rear contact points normalized 
to W, and T indicates the time period required for a dis- 
placement equal to W. From the density snapshots it is 
possible to determine precisely the positions where the 
contact line pins/depins with a resolution that is impos- 
sible to achieve in an experiment. The rear contact line 
gets pinned (snapshots a) before entering the hydropho- 
bic area (yellow region at the base of the drop). A sudden 
jump makes it overcome the local energy barrier and en- 
ter the hydrophilic area (red region at the base of the 
drop) (snapshots b). The front contact line, in turn, 
stops before entering the hydrophobic regions. As the 
drop pins with an increasing advancing angle, it slowly 
penetrates through the hydrophobic area (snapshots c), 
similarly to what happens in Fig. O A sudden jump 
follows when the front enters the hydrophilic area (snap- 
shots d). We have also extended the simulations to the 
heterogeneous surface investigated in [l3| , which presents 
a lower wettability contrast of 10° and found only peri- 
odic oscillations in the velocity without the occurrence of 
a stick phase in agreement with the experiment. 
The bottom graph of Fig. [3] presents the analysis of the 
balance equation (pQ) for a time frame of a period T of 
the stick-slip dynamics and compares it to the homoge- 
neous case, for the same Bo. The gravitational force F g 
is constant in time. When the drop is pinned, F g is al- 
most balanced by F cap (snapshots c). Immediately after, 
the front contact line jumps forward and the drop depins 
(^cap —> 0): the drop experiences a localized acceleration 
a(t) with a consistent dip in the dissipation function. The 
drop gets pinned again until the rear contact line jumps 
forward. The homogeneous case is instead stationary: 
the energy provided by F g is almost entirely transferred 
into dissipation, apart from the deformation of the inter- 
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FIG. 3: (Color online) Top: time evolution in dimensionless 
units (see text for details) of the front (left axis) and rear 
(right axis) contact points of a 2D drop on the heterogeneous 
substrate. The horizontal sequence (on the top) shows the 
whole drop in terms of density snapshots: red (blue) color 
is associated to high (low) density regions. The vertical se- 
quence (on the right) reports the corresponding enlargements 
to better appreciate the position of the contact line. Con- 
sistently, the yellow (red) areas at the drop base represent 
hydrophobic (hydrophilic) regions. Bottom: time evolution 
of the various terms in the balance equation JTJ) for the het- 
erogeneous and the homogeneous case chosen to display the 
same static morphology. 



face which causes a term F cap smaller by a factor ^10 
with respect to the heterogeneous case. Overall, we see 
that the effective dissipation in the heterogeneous case 
is strongly suppressed as compared with the stationary 
homogeneous case. This is because the large wettabil- 
ity contrast causes additional energy to be stored in the 
non-equilibrium configuration of the drop which can pin 
before the contact lines jump forward. 

Our next goal is to quantify the dependence of the drop 
average speed on the inclination angle. Again, this will 
be achieved by comparing the homogeneous and hetero- 
geneous cases. We have performed measurements with 
different coatings of glass slides produced with different 
molecules and methods, obtaining contact angles ranging 
from <9 eq = 71° ± 2° to <9 eq = 115° ± 2°. The sliding of 
water drops of volume V down these surfaces has been 
studied at different inclinations a. For angles above the 
sliding angle, the drop moves at the constant speed U, 
which increases with the inclination angle. Such motion 
is the result of a balance between the in-plane compo- 



4 



nent of the drop weight and the viscous resistance, plus 
a capillary force related to the nonuniformity of the con- 
tact angle along the drop perimeter 0, HjJ . For small 
drops, this leads to the following scaling law between 
the capillary number Ca = 77E//7, where 7 is the sur- 
face tension and 77 the viscosity, and the Bond number 
Bo = (3V/47r) 2 / 3 pg sin a/7, p being the density (30| : 



Ca ~ (Bo - Bo c )/c(0) 



(2) 



where the function c(0) is related to the solution of the 
hydrodynamic flow [31| in the wedge. In principle, the 
angle is the dynamical contact angle. Nevertheless, it 
is plausible to approximate « # eq either when dynamic 
contact angles do not deviate severely from the equilib- 
rium contact angle or when the mean of the advancing 
and receding contact angles is close to eq 0. We 
therefore get c(0) = (1 — cos 2 eq )(0 eq — sin eq cos eq ) 
multiplied by an angular function relating the radius 
of the bottom contact area to the volume [30|. Data 
displayed in the inset of Fig. [4] refer to drops of V w 30 fil 
sliding on both the heterogeneous and the homogeneous 
surfaces with the same static contact angle. The Bo c 
value is deduced by extrapolating the linear fit of the 
sliding points to Ca = 0. It is evident that the two 
surfaces display a different dynamic behavior, the het- 
erogeneous one being characterized by a larger pinning. 
To better compare the sliding on the various surfaces, in 
the top panel of Fig. [H we plot all the data in terms of 
(Bo — Bo c ). In this way, the points lie on straight lines 
passing through the origin, whose slopes depend on the 
surface wettability. The slopes as a function of eq are 
displayed in the bottom panel of Fig. [4| the agreement 
between the angular dependent prefactor of equation 
(|2|) and the experimental data is quite reasonable for all 
the investigated surfaces, including the heterogeneous 
one. This indicates that the effects of the heterogeneous 
patterning can be readsorbed in a renormalized value of 
the Bo Cl representing the increase of the static energetic 
barrier that must be overwhelmed by gravity before 
the drop starts to move, an observation that bears 
similarities with the results discussed in [lij . 
In summary, we have observed both experimentally and 
numerically that a water drop, sliding across alternating 
stripes having a large wettability contrast, undergoes a 
characteristic stick-slip motion. The average speed of 
this nonlinear motion is an order of magnitude smaller 
than that measured on a homogeneous surface having 
the same static contact angle. The slow down is the 
result of the pinning-depinning transition of the contact 
line which causes energy dissipation to be localized in 
time and large part of the driving energy to be stored 
in the periodic deformations of the contact lines when 
crossing the stripes. This suggests that the motion 
of drops can be passively controlled by a suitable 
tailoring of the chemical pattern. Experiments are under 
way to see how the shape of the tiles forming the pe- 
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FIG. 4: (Color online) Top: relation between the Capillary 
number (Ca) and the Bond number (Bo) for both homo- 
geneous (filled symbols) and heterogeneous (open symbols) 
surfaces. All the Bond numbers have been rescaled by sub- 
tracting the critical Bond number Bo c , which is larger for the 
heterogeneous case (see inset). Bottom: the slope obtained 
from the linear fits of the Ca vs. Bo curve is reported as a 
function of the static contact angle. The dashed line is the 
scaling law predicted by J2j) , calculated for small drops sliding 
down homogeneous surfaces with a wedge dissipation as the 
dominant dissipative contribution [30l. \3u\. 



riodic, heterogeneous patterns affects the drop dynamics. 
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